Boundary Slip as a Result of a Prewetting Transition 
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Some fluids exhibit anomalously low friction when flowing against a certain solid wall. To recover 
the viscosity of a bulk fluid, slip at the wall is usually postulated. On a macroscopic level, a large 
slip length can be explained as a formation of a film of gas or phase-separated 'lubricant' with 
lower viscosity between the fluid and the solid wall. Here we justify such an assumption in terms 
of a prewetting transition. In our model the thin-thick film transition together with the viscosity 
contrast gives rise to a large boundary slip. The calculated value of the slip length has a jump at 
the prewetting transition temperature which depends on the strength of the fluid-surface interaction 
(contact angle). Furthermore, the temperature dependence of the slip length is non-monotonous. 

PACS numbers: 47.10.+g, 68.08.-p, 81.40.Pq 



I. INTRODUCTION 

It is accepted in hydrodynamics that the velocity of 
a liquid immediately adjacent to a solid is equal to that 
of the solidi. Such an absence of a jump in the velocity 
of a simple liquid at a surface seems to be a confirmed 
fact in macroscopic experiments. However, it is difficult 
to obtain the same conclusion using microscopic models. 
It has been noticed that, even in case of simple liquids, 
the no-slip boundary condition is not justified on a mi- 
croscopic level. 

Therefore, the conventional condition of continuity of 
the velocity, or the no-slip boundary condition, is not 
an exact law but a statement of what may be expected 
to happen in normal circumstances. While the normal 
component of the liquid velocity must vanish at an im- 
permeable wall for kinematic reasons, the requirement of 
no-slip can be relaxed. In other words, instead of impos- 
ing a zero tangential component of the liquid velocity at 
the solid, it is possible to allow for an amount of slippage, 
described by a slip length b. The slip length for a simple 
shear flow is the distance behind the interface at which 
the liquid velocity extrapolates to zero 



v s = b [d z v(z)] 



wall ' 



(i) 

where v s is the tangential velocity at the wall, and the 
z axis is perpendicular to the surface. The definition 
of b is explained in Fig. 2] It is clear that boundary 
slip is important only when the length-scale over which 
the liquid velocity changes approaches the slip length. 
Therefore, it is not surprising that the slippage effect 
has not been detected in macroscopic experiments. In 
microfluidic devices, however, where the liquid is highly 
confined, the boundary slip is important^. 

Indeed, water flow in capillaries of small diameter with 
smooth hydrophobic walls has been investigated and slip 
at the wall had to be postulated to recover the viscos- 
ity of wateo 3 ^. These results were confirmed by di- 
rectly probing the fluid velocity at a solid surface us- 
ing total internal reflection-fluorescence spectroscopy^ 
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FIG. 1: Definition of the slip length b for a simple shear flow. 



as well as double focus confocal fluorescence cross- 
correlationi. Several indirect methods were also used: 
the quartz-crystal microbalance (QCM)£, the surface 
force apparatus (SFA)2ii2iiiii£, and the atomic force mi- 
croscope (AFM)iii4. The magnitude of the slip length 
b was sometimes greater than lOOnm for partially wet- 
ted wall o 5 ! 10 . In some cases the shear rate did not affect 
the amount of slip in the observed rang o 5 i 14 ; in others a 
strong dependence on the velocity was found 9 - 10 . It was 
also shown that both surface roughness and strength of 
the fluid-surface interactions affect the wall slipiii* 1 ^. 

The no-slip condition can also be violated in more com- 
plex systems. Boundary slip has been suggested for poly- 
mer meltsiiiiiSiii and liquid crystals^ (for the rotational 
motion of molecules). 

The origin of such large slippages remains unclear de- 
spite considerable theoretical effort. On the theoretical 
side, molecular dynamics simulations have shown that 
the molecules can slip directly over the solid due to the 
fact that the strength of attraction between the liquid 
molecules is greater than the competing solid-liquid in- 
teraction 19 ! 20 ! 21 ! 22 ! 23 : 24 . In general, wall slip was found 
on non-wetted surfaces, i.e. when the contact angle is 
large. However, the simulation results were not entirely 
consistent with the experimental data, by predicting a 
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FIG. 2: Slip length b for a binary mixture. 



much lower slip length2Si2& and substantial slippage only 
at large contact angles^. It has also been demonstrated 
that the surface roughness may both reduce^! and in- 
the friction of the fluid past the boundaries. 
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Other ideas invoke the formation of a new phase at the 
wall. The possible source of the surface layer could be a 
gas (lubricant) dissolved in the liquid, forming bubbles 
nucleating at the liquid/solid interfaca^fl*2ii2Si. Experi- 
mental evidence for the formation of a wetting layer un- 
der flow has been found by Tanaka-iiM. The boundary 
layer can have a lower viscosity than the bulk valued. 
Tuning the size and the properties of this layer one can 
obtain large values of the slip length (see Fig. EJ) . 

Indeed, in a sharp interface limit, when the width of 
the interface is much smaller than the width of the slab 
and the thickness of the boundary layer, we can neglect 
the structure of the interface. Then the problem is re- 
duced to the shear flow of two phases (s for "surface" , b 
for "bulk") with viscosities rj s and r] b , respectively, and 
thicknesses 6 and I, respectively (see Fig- - Denoting 
the velocity profiles in the surface layer and the bulk by 
v s (z) and v b (z), respectively, and using stick boundary 
conditions at z — and z = I + 5, we have 



v s (0) = 0, 
v s (S) = Vb(5), 
v b (l + S) = v , 



(2) 



where the second equation is the condition for continuity 
of the flow field at the interface between the phases, and 
vq is the velocity at the top of the slab. Furthermore, we 
use the condition of zero divergence of the shear stress 
tensor at the interface, 



d z (rj(z)d z v(z)) = 0, 



or 



Vsv' s (6) = r) b v' b (S). 
The solution in the bulk 

V b (z) = KV { 1 + - 

\r]s 5 



(3) 
(4) 

(5) 



with k = (rib jris + 1/$) 1 > results in a slip length— 



Vb 



(6) 



According to Eq. JBJ the boundary slip can be observed 
if the viscosity depends on the composition (r/ b r) s ) and 
the less viscous fraction of the liquid wets the walls better 
than the more viscous one (rfc > rj s ). It is also clear that 
there are two ways to obtain a large slip length. First, 
by having a macroscopically thick boundary layer, since 
the slip length has the same order of magnitude as the 
thickness of this layer. Second, by providing large values 
of the viscosity contrast r) b /r} s , e.g. for a gas layer—. 

A more realistic description should allow for a prewet- 
ting transition— i^I for the liquid/gas or liquid/lubricant 
mixture and take into account the structure and the finite 
width of the interface region. 

The aim of this paper is to include these effects and 
relate the slip length to the wettability of the walls, com- 
position of the mixture, and thermodynamic parameters 
of the system. We show that the prewetting transition 
can give rise to a large boundary slip of the fluid by gen- 
crating a thick film of a phase-separated 'lubricant' at 
the wall which has a lower than the bulk fluid viscos- 
ity. Indeed, if we choose typical values, thickness of the 
wetting layer (thick film) 5 w lOnm, viscosity contrast 
Vb '■ = 3 : 1, we find b w 20nm, i.e. the prewet- 
ting film can indeed give a large slip length. This value 
can be further increased. Indeed, if the phase separation 
occurs upon cooling (heating), the thickness of the wet- 
ting layer increases with the increase (decrease) of the 
temperature: for a system with short-range forces in the 
vicinity of the wetting temperature T w , the thickness di- 
verges aaZLS oc - ln(\T w -T\). 

The paper is organized as follows. In Sec. [H] we out- 
line the approach that allows us to calculate the phase 
diagram, order parameter profiles, and slip length. Sec- 
tion II I II gives a summary of the results. Finally, in 
Sec. IIVI we present some brief conclusions. 



II. THEORETICAL MODEL 

Phase separation phenomena in binary and polymer 
mixtures have been intensively studied by theory, exper- 
iment, and simulation. While the most detailed infor- 
mation is available for the static bulk behavior, much 
more interesting phenomena occur when studying the 
dynamics22i2242iii, or the behavior near surfaces and 
in confined geometriea2iii224iMi. The theoretical under- 
standing of phenomena which combine both aspects (i. e. 
dynamics near surfaces) is still at its infancy^, while 
there are many experiments2iSi2&. 

Our theoretical approach splits the problem of shear 
flow of a binary mixture near surfaces into several inde- 
pendent tasks. First, we calculate the equilibrium order 
parameter profiles, completely disregarding the flow. Re- 
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striction to equilibrium thermodynamics allows us to in- 
troduce a change of ensemble: we fix the chemical poten- 
tial difference (semi-grand-canonical ensemble) instead of 
fixing the composition, which is conceptually and com- 
putationally easier. The order parameter profile then re- 
sults in a viscosity profile, which in turn allows us to cal- 
culate the stationary velocity profile by solving Eq. J2J), 
again using stick boundary conditions at both surfaces. It 
is not immediately obvious that this split-up is justified. 
We discuss the restrictions and underlying assumptions 
of this approach in Appendix El Finally, the slip length 
is calculated from the stationary velocity profile. 

A. Free Energy of a Binary Mixture 



which also takes into account a dependence on the overall 
density, which can vary along the profile (see Appendix 
0). 

The term (V</>) 2 is needed to provide spatial struc- 
ture to the theory: at phase coexistence, there are two 
bulk equilibrium order parameter values 4>+ and ef>— with 
the same free energy density, /(</>+) — /(</>-)■ Without 
the gradient term, a structure with very many interfaces 
between the and </>_ phase would be entropically fa- 
vored. The term (fc/2) (V(j>) is the simplest one which 
penalizes interfaces. While this is justified near the crit- 
ical point, where interfaces are very wide and the order 
parameter varies smoothly, a more realistic description at 
strong segregation (where the interface becomes rather 
sharp) would require higher-order gradients, too. 



To describe the bulk phase as well as the interface 
structure, 'phase-field' models 36 are often used. In this 
approach the order parameter <j) is introduced. For a 
binary mixture <f> is a composition variable, defined as 
<f> = (n\ — ni)j(ri\ + 712), where the rii are the number 
densities of the two species. This order parameter varies 
slowly in the bulk regions and rapidly on length scales of 
the interfacial width. The unmixing thermodynamics is 
described via a free energy functional. 

Since the material is confined in a container in any 
experiment, phase separation is always affected by sur- 
face effccts 47 - 4 §i22i£2i£i. To include them, appropriate 
surface terms responsible for the interaction of the liq- 
uid with the container walls are added to the free 
energy 

■3fi.37.42.fi2.fi3.B4 

In the phase-field approach the semi-grand potential of 
a binary mixture is written as^i 



^a 2 (V0) 2 + f(cj>) - ^ 



+ *s, (7) 



where a is a normalization length of the order of the size 
of a molecule, f(<p) is the Helmholtz free energy density 
of the mixture, while /x is the chemical potential thermo- 
dynamically conjugate to the order parameter 4>, and ^ s 
is the surface energy. 

The explicit form of the Helmholtz free energy /(</>) 
varies depending on the type of mixture. However, the 
simple observation that the two phases must coexist im- 
plies that there are two minima in the free energy at 
the respective values of the order parameter. We here 
adopt the mean-field model for a regular (symmetric) 
mixturo 55 ! 56 



kf>T — - — In 



1 + 1 



In- 



, (8) 



where the first term on the right-hand side corresponds 
to the excess energy of mixing. Note that this is one of 
the simplest models to describe unmixing; a more realis- 
tic description would need a more sophisticated function, 



B. Surface Free Energy 

To describe the interaction with the walls we used the 
quadratic approximation for the surface energy22^i 



-heps 



dS, 



(9) 



where <p s is the surface value of the order parameter and 
the parameters h and 7 are referred to as the short-range 
surface field and the surface enhancement, respectively. 
The short-range surface field, ft,, is a measure of the at- 
tractiveness (or repulsiveness, if negative) of the surface 
to the component 1. In real systems it can be of either 
sign and of any magnitude. The surface coupling en- 
hancement, 7, represents the effect that a molecule close 
to the substrate has fewer neighbors than a molecule 
in the bulk; 7 is estimated to be small and negative. 
We have implicitly assumed that all surface effects are 
of short range, thus f s depends on the local concentra- 
tion at the walls only. Because of long range van der 
Waals forces this is not fully realistic. However, for large 
enough wall separations the differences to the short range 
case are rather minor. 

We introduce dimensionless units by setting a = 1, X — 
1, and ks = 1. Hence, energies are given in units of Xi 
temperatures in units of x/^b, and lengths in units of a. 
Furthermore, a value of k — 1 has been used throughout 
the calculations of this paper. 



C. Euler-Lagrange Equations 

In thermal equilibrium the grand potential J7J must 
be minimal. Variation of this functional yields an Euler- 
Lagrange equation together with two boundary condi- 
tions. Due to translational symmetry in x and y direc- 
tion, the problem is one-dimensional. In contrast to the 
situation considered in Introduction, we now focus on 
the case of two identical walls separated by a distance L, 
which is chosen large enough such that the two wetting 
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layers do not overlap, and bulk behavior is established in 
the center of the slab. It is then convenient to choose the 
origin of the coordinate system at the center of the film. 
In this coordinate system the Euler-Lagrange equation 
and boundary conditions read 



d 2 <j> 1 



±k -- + h + 70 

oz 



2 1-, 



• M = 0, 
= 0. 



(10) 

(11) 



z=-L/2,L/2 



This boundary-value problem was solved using the relax- 
ation method 58 . In general, it can have up to three differ- 
ent solutions: one stable, one metastable, and one unsta- 
ble. The relaxation method yields only the metastable 
and the stable solution. The unstable solution with 
the highest free energy has negative response function 
(d(p/dyL)T and is eliminated by the relaxation method 
automatically. 

To select stable solutions, we calculated the grand po- 
tential of the mixture, fl, for both stable and metastable 
solutions and chose the solution with the lowest grand 
potential. This allows the accurate determination of the 
phase diagram. 
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FIG. 3: Typical order parameter profiles. Only the part 
next to the wall is shown. Thickness of the slab L = 200, 
[i — —0.002, 7 = —0.01, h = 0.2. For this chemical potential 
the thin-thick film transition occurs at T w 0.35. 



III. RESULTS AND DISCUSSION 

A. Order Parameter Profiles and Prewetting Phase 
Diagram 



D. Velocity Profiles and Slip Length 

To calculate the velocity profiles, we assumed that the 
viscosity of the system is simply a linear combination of 
that of the individual components 

1 + 6(z) 1 - 6(z) , \ 

ri{z)=Va ^~ L +Vb (12) 

where the viscosity contrast between the two components 
has been chosen as rj s : rjt, = 1 : 3. Of course, a more re- 
alistic description would have to introduce a more com- 
plicated dependence on <f>, and also take into account the 
dependence on the overall density. The stationary veloc- 
ity profile is the solution of Eq. 

, f z dz 
v(z)=v w c / —-, (13) 

where 

r L/2 dz , s 

c= / -pr, (14) 

and we assumed stick boundary conditions at the walls, 
v(—L/2) — and v(+L/2) = v w . The value of the slip 
length was obtained by fitting the bulk region of the ve- 
locity profile (|13|) with a linear regression. 



Typical order parameter profiles at the same value of 
the chemical potential // = —0.002 and a set of tempera- 
tures are shown in Fig. [31 Only the profiles corresponding 
to the stable thermodynamic state are shown. The pro- 
files have a well-developed flat region in the center of the 
film, which indicates that there should not be finite size 
effects for the chosen thickness of the slab (L — 200). 
This flat region is important for defining the slip length, 
see Eq. (QJ. This requires a well-defined linear velocity 
profile in the bulk. 

In order to understand these profiles, one should note 
that the very small value of /i, combined with the low 
temperature, implies that the system is rather close to 
bulk coexistence. Bulk coexistence, however, is charac- 
terized by two equilibrium values of the order param- 
eter with large absolute value and opposite sign. The 
small negative value of /i singles out the negative or- 
der parameter value, while the absolute value is only 
slightly changed. Introducing the surface with preference 
of the other phase, one obtains a profile which is slightly 
bent up near the surface. The prewetting transition (see 
Refs^Si^4) occurs upon increasing the temperature, and 
results in a sudden increase of the surface excess cover- 
age, i. e. it is a first-order transition between a thin-film 
and a thick-film state. This jump is directly observed at a 
temperature near T = 0.35 in the profiles of Fig. EJ Tak- 
ing into account the remarks made in the introduction 
and Eq. ©, we can already anticipate small slip lengths 
for thin films and large slip lengths for thick films (above 
the prewetting transition temperature). Since the thin- 
thick film transition is a first-order transition, one can 
also expect a jump in the slip length at the prewetting 
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FIG. 4: Prewetting phase diagram calculated for h — 0.2, 
7 = —0.01, L — 200. The solid line is the prewetting transi- 
tion line ending at the prewetting critical point. Metastabil- 
ity limits of thick and thin films (spinodals) are shown with 
dashed lines. A, B, C, I, and II are thermodynamic paths 
used in the text. 
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7 = -0.01, L = 200. 
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FIG. 7: Slip length vs. short-range surface field, /i = —0.001. 



B. Slip Length 



FIG. 5: Slip length vs. temperature calculated for several 
values of the chemical potential, h — 0.2, 7 = —0.01, L = 200. 



transition temperature. 

To calculate the thin-thick film transition temperature, 
we followed the metastable solution up to its stability 
limit (spinodal) calculating the grand potential for both 
stable and metastable solutions. The transition line was 
determined from the intersection of grand potentials of 
stable and metastable solutions (thick and thin films). 
Both spinodals as well as the transition line are shown 
in a prewetting transition phase diagram, Fig. The 
solid curve is the prewetting curve ending at the prewet- 
ting critical point and the dashed curves are spinodals or 
metastability limits of the metastable states. 



The temperature dependence of the slip length for sev- 
eral values of the chemical potential (or average volume 
fraction) is presented in Fig. |S] As we anticipated, b(T) 
has a jump at the prewetting transition temperature, 
when the thick film is formed (path I in Fig.0J. Before 
the jump the slip length is small and practically does 
not depend on temperature. After the transition the slip 
length decreases with the temperature increasing, even 
though the prewetting film gets thicker. This is because 
the bulk volume fraction (order parameter (f>) increases, 
giving rise to a decrease in the bulk viscosity (or, equiv- 
alently, viscosity contrast) (see Fig.[3J). 

If we are above the prewetting critical point (path II 
in Fig. 0J there is no jump-like transition, but rather a 
smooth increase of the slip length with the temperature. 

The dependence of the slip length on the chemical po- 
tential is shown in Fig. If the temperature is below T™ 
then the thin film is always stable and the slip length is 
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small (path A in Fig.0J). Once we intersect the prewet- 
ting transition line, the thick film forms with the increase 
of the average fraction of the more wettable phase (val- 
ues of the chemical potential close to zero) and the slip 
length jumps abruptly to higher values (path B in Fig.QJ. 
Above the prewetting critical point T c pw the slip length 
increases monotonically with the increase of the chemical 
potential (see path C in Fig.BJ. 

Finally, the dependence of the slip length on the sur- 
face field h is shown in Fig. [7| It is clear that there is a 
threshold value of the surface field (contact angle) when 
the thick film is formed. At this value of the surface field 
the slip length increases abruptly. Below the threshold 
the slip length is small and practically does not depend 
on the surface field. 



IV. CONCLUSIONS 

To summarize, when the prewetting transition occurs 
in a flow experiment, it may indeed generate a strong slip- 
page. Prewetting provides a mechanism of generating a 
macroscopically thick film at the wall. If this film has a 
lower viscosity than the bulk value, a strong slippage can 
be observed above the prewetting transition temperature. 
The value of the slip length has a jump-like dependence 
on temperature, concentration of the phase-separated liq- 
uid, and surface field (contact angle). 

The mean-field model of wetting is rather general and 
can be applied to liquid-gas systems, binary mixtures, 
as well as incompressible polymer mixtures in the long 
wavelength approximation^. This implies that the large 
boundary slip due to prewetting can be observed in all 
these systems. 

Another goal of this paper is to stimulate accurate 
quantitative measurements of the slip length combined 
with measurements of the fluid wetting properties. This 
will allow to study the underlying microscopic mecha- 
nisms of slippage, and choose an adequate model for ev- 
ery experimental situation. In particular, it would be 
very interesting to measure the slip length as a function 
of temperature while the system undergoes a prewetting 
transition, with the adsorbed species having a lower vis- 
cosity. Our results indicate that in such a case the slip 
length should undergo an abrupt change, and we believe 
that this will be true beyond the limitations of our ide- 
alized model. 



APPENDIX A: HYDRODYNAMICS OF 
SLIPPAGE IN BINARY FLUIDS 

We consider a binary fluid (species 1 and 2) confined 
between two infinite parallel planes to a thin slab of thick- 
ness L. The plane normal is in the z direction. We denote 
the mass densities of the two species by p\ and pi , respec- 
tively, and introduce the linear combinations p = p\ + p2 
(total mass density) and Ap — p\ — p^. 



Concerning the thermodynamics of the system, we no- 
tice that the pressure tensor p a p (Greek letters denote 
Cartesian indices) is anisotropic as a result of the fi- 
nite size effect and the interface contribution to the free 
energy^. Since the system is fluid, there is no elastic re- 
sponse to shear, and hence p a p — for a ^ (3. Further- 
more, pxx — Pyy for symmetry reasons. Apart from this 
anisotropy, we assume that there is no further inherent 
anisotropy. In particular, we assume that all transport 
coefficients (intcrdiffusion coefficient, thermodiffusion co- 
efficient, viscosities, etc.) retain their simple scalar na- 
ture as in the isotropic macroscopic bulk fluid. The vis- 
cous stress tensor is hence written as 

a a(} = j] {dayp + d p v a - ^6 a pd 7 v^j + C<W<9 7 v 7 , (Al) 

where the Einstein summation convention is implied, v 
denotes the fluid flow velocity, while r\ and £ are the 
shear and the bulk viscosity, respectively. We assume 
that these parameters depend on density and composi- 
tion, i. e. r\ = rj(p 7 Ap) and C, = C(p, Ap). 

We consider the situation that the fluid is weakly 
driven such that a shear flow in x direction develops, 
with shear gradient in z direction. In this limit of weak 
driving, it is reasonable to assume that no symmetry is 
broken except translational invariance in z direction. The 
system remains translationally invariant in x and y di- 
rection, and gradients (of any quantity) occur only in 
z direction. We also assume that the system is kept at 
constant temperature throughout, i. e. that the heat pro- 
duction is negligibly small. This is reasonable for small 
shear rates (note that the heat production is proportional 
to the square of the shear rate). In addition, the heat 
conductivity is quite large for many real fluids. 

Under these assumptions, we seek a stationary solu- 
tion of the hydrodynamic equations of motion^ for the 
outlined geometry and symmetry. Since the velocity flow 
field v = (v x , 0,0) is defined as the center-of-mass veloc- 
ity of a volume element, the dynamics of p is governed 
by pure convection: 

d t p + d a {pv a ) = 0. (A2) 

A stationary solution implies dtp = 0, while our geometry 
leads to d a {pv a ) = d z (p(z)v z ) = 0. The continuity equa- 
tion is therefore identically fulfilled, and a non-constant 
mass density profile p(z) is permitted. 

For the density difference there is also interdiffusion: 

d t Ap + d a {Apv a ) = -d a j a , (A3) 

where j is the interdiffusion current. Again, the left hand 
side vanishes identically for our flow. Furthermore, the 
interdiffusion current vanishes in the stationary state: 

3 = 0. (A4) 
The (full nonlinear) Navicr Stokes equation is written as 

9t {pVa) + dp {pV a Vfj) = -dffPaji + dp<r a p. (A5) 
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For our flow, the left hand side vanishes identically, and 
hence 

d z p zz = d z Gzz = 0, (A6) 
= dzp x z(z) = d z a xz = d z (r]d z v x ) , (A7) 
= d zPyz = d z a yz = 0. (A8) 

From Eq. l)A6j) one sees that the pressure profile p zz must 
be constant, while Eq. I|A7(I allows to obtain the velocity 
profile via integration, as soon as the viscosity profile rj{z) 
is known from the profiles p(z) and Ap(z). 

We now turn to the constitutive equation for the inter- 
diffusion current j. In non -equilibrium thermodynamics, 
the dissipative currents are assumed to be linear in the 
gradients of the intensive thermodynamic variables. A bi- 
nary system has three independent thermodynamic vari- 
ables, for which we can take any appropriate set. For our 
purposes it is particularly useful to choose the pressure p, 
the temperature T, and the chemical potential /i, which is 
defined as the variable which is thermodynamically con- 
jugate to the order parameter <f> = (n\ — n2)j(n\ + 712), 
where the rii denote the particle number densities (see 
main text). Macroscopically, the constitutive equation 
would thus read 

j a = -r&n - v 2 d a T - r 3 d aP , (A9) 

where Ti are suitable scalar Onsager coefficients. Note 
that the pressure gradient may appear since p is not the 
variable conjugate to Ap/p (which is often used in the 
literature^), but rather to (j). We now generalize this 
equation to the case of an anisotropic pressure tensor, 
but retain the scalar nature of the Onsager coefficients 
in accordance with our assumptions (similar to what is 
done in model H 61 i 62 ). We thus find 

j a = -Tid a fj, - T 2 d a T - T 3 dpp al3 , (A10) 

or, taking into account that there are only gradients in z 
direction, and that j vanishes, 

= -Txd z p, - T 2 d z T - T 3 d zPzz . (All) 



Now, d z T vanishes due to our assumption of an isother- 
mal system, while d z p zz vanishes as a consequence of 
the Navier-Stokes equation (see Eq. IA6JI . For this rea- 
son, the profile of the chemical potential, p{z), must be 
a constant, too. 

In summary, we find that under the given assumptions 
(isothermal system, isotropic Onsager coefficients, weak 
driving) the conditions for the stationary state are iden- 
tical to those in thermal equilibrium (all intensive vari- 
ables must have constant profiles). This permits to first 
calculate the density profiles just as equilibrium profiles, 
completely disregarding the flow, and then, in a second 
step, to calculate the velocity profile by solving Eq. i|A7(l . 
This has been done in the main text for a simple model 
for the unmixing thermodynamics, and a linear depen- 
dence of the viscosity r\ on the order parameter. Note also 
that the restriction to pure equilibrium thermodynamics 
allows us to introduce a change of ensemble: instead of 
considering the composition as fixed, we rather view p as 
fixed (semi-grand ensemble), which is conceptually and 
computationally easier. Of course, a quantitative com- 
parison with experiments is not possible for such simple 
models; one would have to use much more sophisticated 
free energy functionals, and a better model for the con- 
centration dependence of r\ as well. Furthermore, one 
should expect that at moderate shear rates only hydro- 
dynamic instabilities (e. g. bubble formation near the 
surfaces) should occur which invalidate the assumption 
of translational invariance in x and y direction. 
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